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1. Introduction 



The construction of a quantum theory of fundamental processes that is central to 
physics. Despite this, quantization of even the most elementary systems has proved to 
be far from straightforward. For a system whose configuration space is the whole of the 
real line, R, or any finite dimensional vector space, there are a number of equivalent 
quantization schemes. One approach is the canonical operator formalism in which the 
phase space variables x and p are replaced by operators x and p, and their Poisson bracket 

{x,p} = l (1.1) 

is replaced by the commutator 

[x,p]=i, (1.2) 

the Heisenberg-Weyl algebra. Quantization is completed by finding an irreducible unitary 
representation of this algebra; it may be shown that any such representation is unitarily 
equivalent to the usual one on the space of square-integrable complex- valued functions on 
R, L 2 (R,C) [1]. 

An alternative and very beautiful formulation of quantum mechanics on R is the 
path integral. In this approach, the transition amplitude between an initial state % at time 
U, and a final state / at time tf, is 

G(f,t s ;i, ti ) = / V[ x (t),p(t)]e J: ' ^ l - HM) (1.3) 

In the path integral, V[x(t),p{t)] is the Feynman- Wiener measure on the space of all paths 
in phase space that interpolate between the inital and final configurations. This formula 
applies to any system with Hamiltonian H. and can be shown to be equivalent to the 
operator formulation under certain circumstances [2]. 

In order to treat more complicated problems, we would like to extend the ideas of 
path integration to more general situations such as are encountered in condensed matter 
physics, gravitation and a-models. Our intention is to construct a path integral whose 
perturbation theory respects the global structure of the system. In a previous publication 
[3], we constructed a quantum mechanical path integral for a system whose configuration 
space is R + , the positive real line, and applied it to a minisuperspace model of quantum 
gravity. We showed how the fact that the configuration space variable was restricted to 
be positive led to a different form of the path integral measure than is usually assumed, 
and that the unitarity of time evolution and the question of topology change are crucially 
affected by the correct identification of the path integral measure. 
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It seems clear that related considerations will have a considerable bearing on any 
complete quantum theory of gravity, and this paper is a step in that direction. Specifically, 
we extend the construction of the path integral in [3] to a more general finite-dimensional 
non-compact configuration space of the form G/H, where G and H are two Lie gorups. 
As we shall show, this path integral respects the global structure of the system and allows 
the development of a perturbation theory which goes beyond the conventional weak field 
expansion in analogous field theory systems. 

Although we shall treat quantum mechanics in this paper, our motivation is the field 
theoretic analog to which we intend to return in a future paper. Our primary aim here 
is to develop a perturbation scheme which does not have the drawbacks of the Feynman 
diagram approach usually employed. The first drawback is that conventional perturbation 
theory is based on Gaussian integration, and this means that one allows the integration 
variables to range from minus infinity to infinity which is an approximation that does 
not respect the global structure of general non-linear systems. Secondly, for many non- 
compact configuration spaces, the normal coordinate expansion does not probe the whole 

space: for example the matrix q ^ ^ SX(2,R) cannot be written in the form 

exp(a l T l ) where T % are the generators of SL{2, R). Thus the Vilkovisky-DeWitt approach, 
while reflecting the local geometry, does not allow one to probe the full range of quantum 
fluctuations. 

A further motivation comes from quantum gravity, for which the space of configura- 
tions, positive definite symmetric matrices at each point in space modulo spatial diffeo- 
morphisms, is an example of the field theoretic extension of this work. It is known that the 
weak field perturbation theory for this system is unsatisfactory because it is unrenormaliz- 
able, however the approach adopted here is naturally tailored to a strong field perturbation 
theory [4] . Although this is the most interesting regime of quantum gravity, it is the least 
understood and we expect the work described here to be a step towards allowing us to 
probe this domain. 

We should make it clear that the path integral described here is an alternative to the 
usual configuration space or phase space path integrals which are also correct in the sense 
that they satisfy the relevant Schrodinger equation [2]. We argue, however, that the path 
integral we propose is useful if one wishes to do the sort of "global" perturbation theory we 
are interested in, whereas the configuration space path integral is more adapted to weak 
field perturbation theory . 

The construction we shall present takes as its starting point the canonical quantum 
mechanics of the relevant system, which we shall take to be the finite dimensional coset 
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space G/H where G is non-compact and H is compact. As has been explained elsewhere 
[5], [6], the procedure for finding an algebra of functions on the classical phase space 
closed under Poisson brackets and then turning these functions into operators leads to the 
"canonical" group of the semi-direct product form, WxG : where W is a vector space in 
which G acts linearly with an orbit G/H. 

The space of quantum states corresponds to unitary representations of this group. 
There may be many inequivalent quantizations of this system, in contrast to the case of 
quantum mechanics on R. The inequivalent quantizations of G/H are labelled by irre- 
ducible representations % x of if, with the quantization corresponding to w x realised on 
the Hilbert space of square-integrable functions taking their values in the vector space H x 
carrying the (finite-dimensional) irreducible representation ir x of H, ie. the Hilbert space 
is L 2 (G/H, H x ), [5], [6]. Amongst these quantisations is the "usual" one on L 2 (G/H,C), 
corresponding to the trivial representation of H. 

As a first example, we shall concentrate on the case of n x n positive symmetric 
matrices, Symm + (n,H), that is symmetric matrices whose eigenvalues are all positive 
definite. This space is isomorphic to the coset space GL(n,TC)/0(n) [7] since GL(n, R) 
acts transitively on Symm + (n, R) by conjugation, 

Y^g T Yg Y e Symm + (n, TL),g e GL(n, R) , (1.4) 

with little group at the identity, 0(n). The canonical group [7] is Symm{n, R)xGL(n, R) 
with Lie algebra relations 

[Yab, Y c d] = 

[n a \n c d ]=i(5 d a n c b -5 b c n a d ) (1.5) 
[Y a b,7r c d } =i(6iY bc + 8iY ac ), 

where Symm(n, R) is space of all real n x n symmetric matrices, and is isomorphic to 
R2 n ( n + 1 ). The natural GL(n, R) invariant Hamiltonian on Symm + (n, R) is 

£ = \^a b ^b a ■ (1.6) 

In this paper, we shall treat only the quantisation in which 7i x = C and in this case the 
representation of n and Y is given by 

7, a ^{Y) = -iY ab ( iP(Y) 

\ dY Jbc . (1.7) 

Y ab i;(Y) = Y ab VCO 
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with 

= (1 + <W(tJ-) (no sum) (1.8). 

afe \OYabJ 

The plan of this paper is as follows: in section 2 we outline how to carry out Fourier 
analysis on Symm + (n, R) which we then use to construct a formal path integral. In section 
3 we shall show how this form of the path integral, rather than the configuration space or 
phase space form, leads naturally to the "global" perturbation theory; we then calculate 
the first order in perturbation theory for the quantum mechanics of a particle moving 
on Symm + {n, R) in a simple potential, in order to demonstrate that the terms in the 
perturbation expansion may be evaluated in detail. The form of this expansion suggests 
that this highly nonlinear model has an exact solution. We verify that this is indeed the 
case. In section four, we show how this construction can be generalized to more general 
non-compact symmetric spaces. Finally, some conclusions are presented in section 5. 

2. Fourier-Helgason Transform and Path Integral 

The Fourier modes e ipx on the real line R have two important properties. The first is 
that they are "complete" in the sense that the Fourier transform f(p) of a function f{x), 
given by 

/oo 
f{x) e ipx dx, (2.1) 
-oo 

may be inverted by the Fourier inversion formula 

f(x) = ^ fip) e~ tpx dp, (2.2) 

and secondly that they are eigenfunctions 1 of the invariant differential operators on R, 
the most important example of which is the Laplacian -jjjp. 

These two properties are extremely useful since they allow one to construct the prop- 
agator for a free particle, that is the propagator for the Schrodinger equation with free 
Hamiltonian — It is 

/oo 
dpe ip(x f - Xi ) e kiv\t f - ti ) (2 _ 3) 
-oo 

1 Since e ipx is not normalizable with respect to the usual L 2 inner product, it is not, 
strictly speaking an eigenfunction; this problem may, however, be treated by rigged Hilbert 
space techniques [8]. 
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These Fourier modes are also used in the construction of the path integral for a particle 
moving on R. In order to find the phase space path integral expression for the propagator 
(xf,tf | Xi,ti) one inserts complete sets of states | x) and | p) at N intermediate intervals. 
| x) is an eigenstate of the operator x with eigenvalue x and | p) is an eigenstate of p with 
eigenvalue p. They satisfy the completeness relations 

/oo 1 Z" 00 

dx\x)(x\=— / dp\p){p\=l. (2.4) 



Thus 



Xf - X N ) 



(X f, t f \Xi, U) ~ \im^ ^JJ J dx r(^^j J dP r^j <K 

expz j ^2 (^Pr+i(x r+1 - x r ) - At H(x r , P r +i)^ j 



so that, were the limit to exist, 



(2.5) 



{Xf,tf\Xi,U) 

N r-OO / -J^ \ N r-OO \ 

Y\ J dXr \2TxJ J dPr J { x f,tf\ x N,tN}{xN,tN\PN,tN}--- 
. . . |a? r -|_i , £ r -|_i) (x r _|_i , t r -\-\ \ P. r-\-\ , t r ^i) (P. r-\-\ i ? ^r) (^n | • • • 

. . . (P 1 ,t 1 \x i ,t i ). 
We may now use the fact that 

and 

(P r+1 ,t r+1 \x r ,t r ) = e - iP r+lXr-i(tr+l-t r )H(x r ,P r+1 ) 

to order (t r+ \ — t r ) and hence implicitly to lowest order in %. 
Therefore 



(2.8) 



(x f ,tf\xi,U) ' = ' Z j VP j Pxexpji^ / dt(P(t)x(t) -H(x,P))^j (2.9) 



where Z is a normalisation constant. 
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This last equation is the usual starting point for perturbative calculations. If H = 
Hq + XV with some Hamiltonian H Q for which the path integral is known, and potential XV, 
then by standard methods one can calculate the path integral perturbatively as a series in 
positive powers of A. However, since the limit N — > oo does not really exist, one can argue 
that the path integral is largely a mnemonic and is correct in as much as it reproduces this 
expansion correctly. This is not to say that the path integral is devoid of content; indeed 
the path integral can be used to treat various phenomena non-perturbatively, an example 
is provided by the use of instantons to treat tunneling processes. The expansion of the 
propagator Gy for H = H + XV in positive powers of A is most easily derived from its 
integral equation [9], namely 

r-tf r-OO 

G v {xf,tf,Xi,U) = G (xf,tf,Xi,ti)-X dt dx G (xf,tf,x,t)V(x)G v (x,t;x i ,t l ). 

Jti J — oo 

(2.10) 

Thus the perturbation series may be written, schematically, as 

G V = G + XG VGo + X 2 G VG VG + (2.11) 

where 

ftf r°° 

G VG = dt dxG (x f ,t f ;x,t)V(x)G (x,t;x i ,t i ) , (2.12) 

Jti J — oo 

and so forth for higher order terms. 

We have set out this well-known derivation in detail since it may be repeated in much 
more general situations such as the general (non-trivially induced) quantisations of motion 
on the arbitrary finite-dimensional coset spaces G/H in which we are interested. 

Rather than treating the general case here, we shall start with the example of motion 
on the space of n x n positive symmetric matrices Symm + {n, R). The beginning of this 
section follows Terras [10] closely. The analogs of the Fourier modes in this case are the 
so-called power functions p s (Y[k]), where Y E Symm + (n, R). These functions are labelled 
by s G C n , k E 0(n) and Y[k] means conjugated Y by k, thus Y[k] = k T Yk. Then 

n 

p B (Y) = Y[\detY j \ a * , (2.13) 

where Yj is the j x j matrix obtained by taking the first j rows and columns of Y. Power 
functions are eigenf unctions of the invariant differential operators L. For GL(n, R), a basis 
for the commutative ring of invariant differential operators is given by 

L (r) = Tr 7r r (2.14) 
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with r = 1 . . . n. For any invariant differential operator L, the eigenvalues depend only 
on s 

L p s (Y[k]) = X L (s) p s (Y[k]) . (2.15) 
The Fourier- Helgason transform /(s, k), of a function f(Y) on Symm + (n, R) is given 

by 

/(s, k)= f duty) f(Y) p~^Y[k\) (2.16) 

where d[i(Y) is the measure on S't/mm + (n,R) invariant under GL(n, R), and the bar 
indicates complex conjugation. The inversion formula is 

f(Y)=u n [ ds f dk f(s,k) Ps (Y[k}) \c n (s)\- 2 (2.17) 

Js=v+iu JkeK/M 

where 

and each component of u is real and runs from — oo to oo. 

= ft ^rrn < 2 - 19 > 

and the Harish-Chandra c-function is given in terms of the beta-function 

*•■»>- ^8 ■ ^ 

as 

-w= n B(i ";Tr ( ; j+ < l ( i ) r +1)) - < 2 - 2i > 

l<i<j<n-l D \2i2\J L ^ L )> 

Finally, the space over which the fc's are integrated is the coset space K/M where K is 
0(n) and M is the set of diagonal matrices with non-zero entries of ±1. That this integral 
is over K rather than K arises because of the degeneracy in specifying k given a particular 
power function coming from a specified Y. In other words M is the automorphism group 
of the invariant power functions. By convention, one normalizes the k integration over 
K/M to be unity 

dk = l (2.22) 

K/M 
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For simplicity we shall write the above inverse transform as 



f(Y) = J d» He ig(a,k) f(s,k) Ps (Y[k)) (2.23) 

As in the case of quantum mechanics on R, we may use this transform to give an expression 
for the free propagator where the Hamiltonian is the invariant Laplacian on Symm + (n, R). 
The propagator Go(Yf, tf; Y t , U) is the solution to 

H G (Y f ,t f ; Y u U) = i—G (Y f ,tf, Y,U), (2.24) 

the Schrodinger equation, with H = ^Tv a b ^b a and lim ti _ >t/ Go(Yf, tf, Y i: ti) = 8(Yf,Yi), 
and it is given by 

G (Y f ,tfY i ,t i ) = J dfi Helg (s, k)p s (Yi[k])p s (Yf[k]) e ^«o (-)(*/-**). ( 2 .25) 

We may also use the power functions to form a path integral for this system just as in the 
case of quantum mechanics on R. We insert complete sets of states, 

J dfx inv (Y) \Y)(Y\ = J dfi Helg (s,k) \s,k)(s,k\ = 1 (2.26) 

to produce the following expression for the propagator associated with the Hamiltonian 
Ho + XV 

(Y f ,t f \Y i ,t i )=[ ( TT d» Hel9 (s n ,k n )d» inv (Y n )) eEl.^^^W^^-W)) 

J Vn=i / (2.27) 

We note that this is not the usual phase-space path integral. However, expanding in pow- 
ers of A gives the usual canonical perturbation expansion in which all integrals are over 
measures and ranges appropriate to the global structure of the problem under considera- 
tion. 

3. The Perturbation Expansion 

We have argued that if we want to develop a perturbation theory that respects the 
global structure of the space of configurations, we should use (2.27) rather than the usual 
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phase space path integral, which is tailored to perturbation theory using Gaussian integrals 
and hence only respects the local geometry. The usual expression does, at least, have the 
benefit of allowing calculation of the perturbation series. We now show that the terms in 
the perturbation series based on (2.27) may also be calculated even though at first sight 
the relevant integrals look intractable. 

For definiteness we consider the case of 2 x 2 positive symmetric matrices and the 
Hamiltonian H + Alndet(Y), where H is the free Hamiltonian, given by the invariant 
Laplacian. We now evaluate Gy to the lowest order in perturbation theory in A. 

The first order term in the expansion is thus 

A^ ' dtj d/j inv (Y) Go(Y f ,t f ;Y,t) lndet(Y) G {Y,t;Y h U) . (3.1) 
We show that this is equal to 



-(t / -t,)(lndet(Y l )+lndet(Y / )) ) Go(Y f ,t f ;Y it U) 



(3.2) 



Thus with 
it follows that 

U ~ b ~ a 

-"0 = o^a 7Tb 



Y = 



Y n Y 

Y\2 Y 22 
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(3.3) 



2Y°+Y 12 ° 



dY u 



dY u 



+ [2Y 12 -— + Y 22 ' 



dY u 



d d 
2Y 12 —-+Y 22 - 



dY u 



'dY u 



'dYt 



12 



y u » +2 y 12 8 



dY] 



12 



d d 
r \2Y 22 —-+Y 12 



dY 22 

2 



dY 



22 



<9Yi 
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Thus since the eigenfunctions of the Laplacian are 

p s (Y) = Y 1 s ?(Y ll Y 22 -Y 1 2 2 y> 
it follows that the spectrum of the Hamiltonian specified by 



(3.4) 
(3.5) 



H 0Ps (Y) = X Ho (s)p s (Y) 



(3.6) 



is 



Ah (s) = -2s? - 4sl - 4s!S 2 - si 



(3.7) 
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Thus the propagator is given by 

Go(y,t;Yi,ti) = JdsJ dk\c(s)\- 2 Ps (Y[k]) Ps (Y 1 [k]) e -*W«0(*-*i) (3 8) 

The Harish- Chandra c-function evaluated along the relevant contour is 

|c( h iui, — + iv,2)\~ 2 = ttwi tanh7r-ui . (3-9) 

Now, given the form of d^Heig 

= J dfi Helg A [p B (y[k])p t (Y 1 [k]) e -^W(*-*i)^ 
= -i y e^He* 5 ^lndet(yi) - lndet(y) - 4i(* - ^(si + 2s 2 )^) (3 - 10) 

xp 8 (y[fc])p s (yi[fc]) e -^(s)(t-t!) _ 

But along the contour used in the inverse Helgason transform, 

«i^ = (-85^ + 5*^. (3-") 



and 



d 



12 



+ 1 - 2 5hTd^) )p -« |fcl) - 



hence 



(3.12) 



lndet(y) + 4i(t - ^)((y) n — — - + (y,) 22 ^^+ 

d(y)n o(yj22 
(^)i2^|yH) G (y,t;y,t,) • 

(3.13) 

This allows us to do the Y integration using the fact that 

J dfi inv (Y) Go(Y f ,t f ;Y,t)G (Y,t;Y i ,t i ) = Go(Y f ,t f ;Y i ,t i ) . (3.14) 

11 



Hence we find that 



G 1 (Y f ,t f ;Y l ,t l ) = \ T dt\(\ndet(Y i ) + 4i(t-t i )((Y i ) 11 —^— + (y i ) 227 ^ 7 ^_ 



d(r,)22 

G {Y f ,t f ,Y h U) , 

(3.15) 



'd(Y t ) 12 
leading to 

G 1 (Y f ,t f -,Y i ,t i ) = \(^(t f -t i )(lndet(Y^ + lndet(Y f ))^ G^Y^tf^U) , (3.16) 
as in eqn (3.2). 

It may be seen that the techniques used in this example allow the calculation of the 
perturbation expansion for a large class of interactions, leading to expressions of the form 

G n {Yf,tf-Y l ,t i )=\ n f n {Y f ,tf-Y l ,U) G^Yf.tf-Y^U) . (3.17) 

where f n is a multiplicative function. 

As a last point in this section, we remark that, our expression for the first order term 
for the perturbation expansion for motion on positive symmetric 2x2 matrices in the 
potential V(Y) = Alndet(F) suggests the following Ansatz for the exact propagator for 
this system: 

G V (Y f , tf, Yt, U) = G (Yf, tf, Yi, t j ) e (^(*/-^)(lndet(l-)+lndet(y / ))-^A 2 (t / -t0 3 ) . (gjg) 

That this is the correct exact expression may be checked by direct substitution into the 
relevant Schrodinger equation. We have thus added to the relatively small number of 
quantum systems whose propagator may be found in closed form. 2 

4. General Non-Compact G/H. 

The same methods that were applied to quantum mechanics on positive symmetric 
matrices can be applied directly to the general case of non-compact G/H, at least for the 
trivially induced representation of H , and provided that if is a maximal compact subgroup 
of G. The ingredients of this construction are a knowledge of the invariant Laplacian L^ 2 \ 



2 We believe that there is an error in Feynman and Hibbs [9] p64 equation (3-62) and 

/?2ry-i3 j?it-i3 

that the correct last term in that equation should be —4^- rather than —^24" 
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its eigenfunctions, the Helgason transform and its inverse. These building blocks may be 
found in explicit detail in Helgason [11], Barut and Ragzka [12] and Terras [10]. The result 
is essentially the same as equation 2.27, namely 



This is the transition amplitude between the initial state described by an element of G/H 
given by Fj, and a final state given by Yf. As before, Y[k] is Y e G/H conjugated by 
k E H, and p s (Y) are the power functions associated with the Laplacian L^ 2 \ which is 
taken to be a multiple of the free Hamiltonian H with spectrum A# . In the path integral 
itself, d[ii nv is the canonical invariant measure on G/H, that is the volume form derived 
from the same metric that determines L^ 2 K d^Heig is the measure for the inverse Helgason 
transform. This expression is quite complicated, but Gindikin and Karpelevic [13] have 
given an explicit form for the measure for such G/H. We do not record the specific forms 
here, but refer to original literature on this topic. 

5. Conclusion 

We have constructed a path integral for quantum mechanics on the space of positive 
symmetric matrices whose natural perturbation theory respects the global structure of 
the system. A more usual approach to the problem would be to use the phase-space or 
configuration space path integral. In the former case, one has (schematically) 



and it can be shown [2] that this expression satisfies the Schrodinger equation, so it must 
be equal to our expression. However, although the above equation is correct it is not clear 
how to use it to develop the sort of expression we are interested in. 

An alternative way to treat many non-linear systems using a phase-space path inte- 
gral is to follow Faddeev [14] and embed such a system in a vector space (for which our 
expression becomes the usual one) and use the constraints to restrict integration to the 
submanifold of interest. Unfortunately, this avenue is not open to us in the class of exam- 
ples we have treated here since 5|/mm + (n,R), rather than being defined by an equality 
(eg. a 2-sphere is defined by x 2 + y 2 + z 2 = 1 in R 3 ) is defined by inequalities (eg. the one 
dimensional case is just the positive real line x > 0) . 




(4.1) 



e l E! =1 ( A «o( s »)+ A1/ ( y »)) At 




(5.1) 
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A further way one might have tried to treat the system dealt with in section three 
would be to realise that every positive symmetric matrix is the exponential of an arbitrary 
symmetric matrix, and that the exponential map is one-to-one and onto, so that we could 
have equivalently treated the space of arbitrary symmetric matrices which is a vector space. 
Apart from the fact that in a field theoretic version of this problem such a transformation 
would be extremely delicate (even e^ x ^ is ill-defined for usual scalar field theory), this 
isomorphism is not helpful in the case at hand. The reason is that, although the space of 
nxn symmetric matrices is isomorphic to Ri n ( n+1 ) topologically, the natural metric on the 
former (which is the one we are using) has constant negative curvature. Furthermore the 
Hamiltonian on Rl n ( n+1 ) in which we are interested is the one invariant under GL(n, R) 
and so rather than the usual Laplacian — on R§ n ( n+1 ), we need a quantisation 

adapted to this Hamiltonian. In other words, in order to treat the Hamiltonian of interest, 
we have to use a "strange" quantisation of R5 n ( n+1 ) using Symm(n, H)xGL(n, R) rather 
than the Heisenberg-Weyl group. In fact, viewed as a theory on Rl n ( n + 1 ) j our system 
for the free Hamiltonian is related to one of the integrable quantum systems discussed for 
example by Olshanetsky and Perelomov[15]. 

Lastly we would like to put forward a path integral for a field theory taking its values 
in positive symmetric matrices. By analogy with the quantum mechanical case, we suggest 
that the relevant expression should be 



e t J dx(\ L (s(x))+V(Y(x))) 

and that perturbation theory based on this expression will respect the global structure of 
the configuration space. 

An analog of formula 5.2 is an alternative solution to the problem of finding the 
phase space path integral for quantum gravity. This quiestion has been addressed in a 
large number of publications and often leads to measures of the form d/iLiouvUie and hence 
to configuration space measures dfXi nvar iant- It is then often argued that in dimensional 
regularization the form of the measure is unimportant and hence that finding the correct 
form of the measure is fruitless. As was clear from the example in section three, this is not 
the case for the path integral put forward here. 



G=H[ (e^ 5 (s(aO,/^ 




(5.2) 
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